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This is a model study for the emergence of superconductivity in ferromagnetically ordered phases 
of cubic materials whose crystal structure lacks inversion symmetry. A Ginzburg-Landau-type 
theory is used to find the ferromagnetic state and to determine the coupling of magnetic order to 
superconductivity. It is found that noncentrosymmetricity evokes a helical magnetic phase. If the 
wavelength of the magnetic order is long enough, it gives rise to modulations of the order parameter 
of superconductivity, both in modulus and complex phase. At magnetic domain walls the nucleation 
of superconductivity is found to be suppressed as compared to the interior of ferromagnetic domains. 



I. INTRODUCTION 

If the unit cell of a three-dimensional crystal is non- 
centrosymmetric, i.e., not invariant under the parity 
operation, the spatial inversion symmetry is broken. 
Noncentrosymmetricity allows for the Dzyaloshinskii- 
Moriya interaction that entails the breaking of spin ro- 
tation symmetry due to antisymmetric spin-orbit cou- 
pling (ASOC). Many noncentrosymmetric material fea- 
ture unusual properties, when ordering phenomena such 
as superconductivity or magnetic order, that break other 
symmetries in addition to inversion and spin rotation, are 
present. For example, ASOC and magnetic order lead to 
a helical magnetic structure in MnSi^ and Ba2CuGe207'2l, 
where they give rise to magnetic field-induced ferroelec- 
tricitjS'. 

Superconductivity with a noncentrosymmetric crystal 
structure generically appears in a mixed parity state^. 
In addition, many of the know n no ncentrosymmetric su- 
perconductors, such as CePtaSPSi, CeRhSiiP, and UlP, 
show states of magnetic order and in some regions of 
the phase diagram superconductivity coexists with the 
magnetic order. A simple consideration of the combined 
action of ASOC and time-reversal symmetry breaking 
magnetic fields or magnetization on the energy bands re- 
veals the possibility of spatial modulations of the order 
parameter of superconductivity (finite-g-pairing)2H121, 

The aim of this paper is to study noncentrosymmet- 
ric systems which show coexistence of magnetic order 
and superconductivity. To be concrete, we will focus 
on a cubic crystal without inversion center described by 
the point group O and restrict our study to the case 
of ferromagnetic order, where the wavelength of mod- 
ulations of the magnetization is much longer than the 
lattice spacing. We will use the generalized Ginzburg- 
Landau-approach. In sec. II we determine the nature of 
the helical ferromagnetic phase in presence of ASOC. In 
sec. Ill we subsequently study, how the superconducting 
state is altered on the background of this magnetic order. 
We find that noncentrosymmetricity causes the magnetic 



moment to follow a helical modulation and introduces a 
new length scale for the superconducting order parame- 
ter. Depending on the ratio of the magnetic wavelength 
and this length scale, the order parameter of supercon- 
ductivity either remains homogeneous or exhibits a mod- 
ulation both in complex phase and absolute value in this 
magnetic phase. Finally, in sec. IV we consider a limit 
where the magnetic state resembles a filamentary domain 
structure. We show that superconductivity will nucleate 
in the interior of the domains rather than at the domain 
walls. 



II. MAGNETIC STATE 

In order to determine the magnetic state, the free en- 
ergy density is expanded in the magnetization m{r) as a 
three-dimensional order parameter. The expression must 
be invariant under time reversal symmetry and under the 
symmetry transformations of the cubic point group O. 
Spin-orbit coupling ties rotations of spin and orbital de- 
grees of freedom together, such that the free energy has 
to be invariant under a simultaneous rotation in both 
spin and orbital space. The magnetization m{r) belongs 
to the irreducible (vector) representation r4 of the point 
group O. The second-order terms of the free energy read 



f!^^ := / d?'r\ am? 



+ T3 



+ 



■{■ 



E 



To +Tl 



To - Tl 2 

— 7, — (V X m) 



T2 (V • m) 



\i5ij{dimi)'^ - {dimi){d.jmj)\ 



(1) 






dm ■ (V X m) \, 



where d and r^ (i = . . . 3) are temperature-independent 
phenomenological parameters and a has the usual tem- 
perature dependence a — a'{T — Ty^), with a' > 0. Here, 
T-^ is the transition temperature for a state of homoge- 
neous magnetization. The effect of noncentrosymmetric- 
ity manifests itself in the presence of a Lifschitz- invariant 



term proportional to ■(^12121, which would not be allowed 
had we considered the centrosymmetric cubic point group 

(2) 

Oh- We demand that F^' is bound from below towards 
arbitrarily strong modulations of m(r). Amongst others 
this leads to the condition tq > 0. 

(2) 

The vector field m{r) that minimizes F^^ is given by a 
helix with a wavevector fco- Depending on the parameters 
Ti (i = . . . 3), fco is either ahgned parallel to a coordinate 
axis (6-fold degenerate) or parallel to a body diagonal (8- 
fold degenerate). Let us assume feo = ^o^z from here on. 
In that case, the magnetization takes the form 



m{r) = niQ [sin{zLkoz)ex + cos (fcpz) By] . 



(2a) 



The two signs correspond to two energetically degenerate 
chiralities and fco = \'d\/{2To). The magnitude mo of 
the magnetization would be determined by terms of the 
order m*. For the time being we ignore these terms and 
consider mg as a parameter. The transition temperature 
to the magnetic state is given by 
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and is always larger than Tj^, owing to tq > 0. In the 
limit of a centrosymmetric system (i? — > 0) we recover 
a state of homogeneous magnetization with the corre- 
sponding transition temperature, Ty^. 



III. SUPERCONDUCTIVITY IN THE 
MAGNETIC PHASE 

In this section, we study the influence of homogeneous 
and helical magnetic order as given by eq. (2a) on the 
superconducting order. We assume the temperature to 
be sufhciently below the transition temperature to the 
magnetic phase, such that an emergent superconducting 
order parameter does not change the magnetization con- 
siderably. For simplicity, a complex scalar superconduct- 
ing order parameter ri{r) belonging to the Fi representa- 
tion of the point group O is considered. The expansion 
of the free energy up to the order -q"^ reads 
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where Z? = V — 2ieA/{hc) with A being the vector po- 
tential. Here, b, u and v are temperature-independent 
phenomenological parameters and a has linear temper- 
ature dependence a — a'{T — Tg^) with a' > 0. T^q 
denotes the critical temperature of superconductivity for 
m — 0. Demanding that F^q is bound form below if 77 
has strong spatial fluctuations necessitates 6 > 0. The 
term proportional to u represents the paramagnetic de- 
pairing effects of the magnetization on the Cooper-pair 
formation, thus we assume u > 0. Noncentrosymmetric- 
ity is again reflected by the presence of a Lifschitz-term 



proportional to v, which is forbidden in case with inver- 
sion symmetry. Note that the Lifschitz-term introduces 
a new length-scale |^l| for the superconducting order, 
where ^l '■= b/lvmo) and is not singular at the transi- 
tion to the superconducting phase. In the following we 
shall assume the limit of strong type-II superconductiv- 
ity in the sense that the length scales fc^^^, ^l, and the 
coherence length ^/b/a are assumed to be much smaller 
than the magnetic penetration depth XU^ In this limit, 
we can neglect the effect of the vector potential A and 
replace Z> — >■ V in eq. ([s]). 



A. Homogeneous Magnetization 

Before addressing the helical magnetic order, we first 
consider the effect of a homogeneous magnetization m = 
moBy on the onset of superconductivity, i.e., we take the 
limit fc^ — >■ 00. Minimizing the functional ([3]) with re- 
spect to 77 straightforwardly yields a modulation of the 
order parameter as 



7;(x) = Tyoe*?'/^^ 
The transition temperature is then given by 
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and shows the advantage of the phase modulation of the 
order parameter as compared to the transition temper- 
ature to a homogeneous superconducting state Tg^™ = 
Tqq — urn^la' . The appearance of these phase modula- 
tions of 77 in noncentrosymmetric superconductors in a 
homogeneous magneti c fiel d was already pointed out in 
various other studies^HISIIll, 



B. Helical Magnetization 



We now turn to the more subtle effect of the helical 



magnetization given by eq. (2a) on the superconducting 



state. The variation of F^q with respect to 77 yields 

= [dl + A(B) + 2B cos(2C)] VkAO, (5a) 

where a Fourier transformation of the superconducting 
order parameter in the coordinates x and y was per- 
formed and {kx,ky) = (fccos(/3, fcsintp). 
We introduce here the parameters 
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(5b) 
(5c) 
is substituted by ^ as C •= 



and the z coordinate 
(fcoz - ip -7r/2) /2. 

In eq. (5a) we identify Mathieu's differential equa- 
tion with the variable C- Mathieu's equation cannot be 





FIG. 2. (Color online) Schematic picture of the superconduc- 
tivity order parameter rj^ (r) given by eq. ([7| that nucleates in 
the presence of the sketched helical magnetization. It exhibits 
an amplitude modulation in the z-direction, and in addition 
to that features a phase modulation perpendicular to the z- 
direction (depicted by the colour gradient). 



FIG. 1. (Color online) Portion of the stability diagram of 
Mathieu's equation. Equation ( |5b[ ) defines a parabola At{B). 
As the temperature decreases, the parabola migrates from 
left to right. The transition temperature is reached when the 
parabola touches the region of stability (dots on dashed curve) 
at some critical value Be, which in turn determines k. If the 
curvature of the parabola at fc = is larger than that of the 
stability region (dotted curve), there will be no contact with 
the region of stability until the transition to the homogeneous 
state Be = 0. 



solved analytically in a closed form. Hovifever, one can 
identify a region in the the A-S-parameter space where 
the solutions of cq. (5a) are bounded for all ^ S M. 
This stability region is displayed in fig. IT] as a shaded 
area. For parameter values {A, B) inside this stability 
region, the system would be in a stable superconduct- 
ing state. Equation (5b) defines a parabola At{B) of 



energetically degenerate parameter values in the A-B- 
plane. These parabolas can be labeled by temperature 
T via the temperature dependence of A. For sufficiently 
high temperatures these parabolas do not intersect with 
the region of stability. As the temperature is lowered, 
the superconducting instability occurs when the parabola 
At{B) touches first the boundary of the (shaded) sta- 
bility region in figjl] Different types of solutions de- 
pend on the curvature of the parabola around B — Q 
which should be compared with that of the boundary to 
the stability region, which can expanded at _B = to 
A{B) « -B^/2 + 7B^/12^. If the curvature of At{B) 
is larger than that of the stability region, corresponding 
to the condition 



laifco > V2, 



(6) 



then the touching point is at Sc = (see dotted line 
in Fig. IT]) such that with fc = the superconducting 
order parameter is homogeneous. For the condition op- 
posite to eq. ^ the touching point is at finite values of 
Be « ±{16(2 - ilkD/iy/'^ (< 1) yielding a modulated 
order parameter with a finite fc, as shown by the dashed 
line in fig. [I] Equation ^ can be seen as an analogue to 
the condition on the Ginzburg-Landau parameter k, that 



appears in the discussion of the vortex phase of supercon- 
ductors. As in our case, we compare a superconducting 
length scale (^l) with a magnetic length scale (fc(J^) and 
obtain an inhomogeneous superconducting state if the 
superconducting length scale is shorter. 

For the inhomogeneous case, the solution features a 
continuous circular degeneracy in the fc^^-fcy-plane, pa- 
rameterized by Lp. The elementary solution is given in 
terms of the lowest order even Mathieu function ceo^ by 



'7V (^) = % exp [ifcoCL^c {x cosLp + y sin ip)] 
X ccq [k^z — If — 7r/2, B^) . 
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This form of the order parameter is sketched in fig. [2] 
together with the helical magnetic state. The general 
solution is a superposition yy^ :— 'Y^^=i^i''lipi with the 
complex- valued coefficients C; . To find the set of coeffi- 
cients ci and phases Lpi that mininrizes the free energy, 
we use Abrikosov's parameter /3 defined by 



/3:= 






(8) 



where (. . .) is the spatial average. The solution that min- 
imizes j3 is realized. A minimum requires that the par- 
tial derivatives 9i^,/3 and dciP vanish for all I. For the 
phases, this yields the condition ipi — ipj = nNij/2, yi,j 
with some Nij G 1. An explicit evaluation of /3 for the 
remaining cases n = 1 
to be 



. 4 reveals the optimal solution 



'7^'"*''^ = r\^ + inv±Txl2- 



(9) 



From this result we see that the order parameter ac- 
quires a long-wavelength phase and amplitude modula- 
tion perpendicular to the wavevector of the helical mag- 
netic order feg and an amplitude modulation with the 
same wavevector as the magnetic order. To lowest order 
in B^ (X {2 — fco^L)' the transition temperature of this 
inhomogeneous superconducting state is given by 
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and the superconducting order parameter has the ap- 
proximate form 
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(10b) 



IV. MAGNETIC SOLITONS AND 
SUPERCONDUCTIVITY 

(2) 

The analysis of the free energy F^ ' showed that ASOC 
favors an inhomogeneous magnetic state. So far we ig- 
nored in the discussion the exphcit form of the fourth 
order terms in the free energy expansion. Besides fix- 
ing the magnitude of the magnetic moment, these terms 
may also introduce additional features such as crystal 
anisotropy of the magnetic moments. The fourth order 
terms allowed within the point group O, ignoring any 
gradient terms, read 
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Here, /3o and /3i are temperature-independent phe- 
nomenological parameters. The expression shows that 
not all directions of the magnetization are degenerate in 
energy. If we simply insert the helical solution ( 2a ) we 
obtain 
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which does not minimize this part of the free energy. It 
is instructive to solve for the magnetization m{r), that 
minimizes the total magnetic free energy F^' + F^' , 



in order to understand the effect of the anisotropy term 
/3i qualitatively. As before, we assume that the mag- 
netization has no component in Cz-direction and will 
be a function of z only. With the complex notation 
M{z) = mx{z) + imy{z), the variational equation reads 



= {-Tgdl + idd^ -f a) M{z) 

+ 3(/3o + Pi) \M{z)\^ M{z) + /?iM3(z). 
To linear order in /3i, the solution is given by 
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where now mg = a'iT^^ - T)/[3(/3o + /3i)]. Thus, the 
anisotropy term pins the magnetization m(r) parallel to 
Bx or By (parallel to e^ ± By) for /3i < (for /3i > 0), be- 
sides a modulation of the amplitude |»Ti(r)| — \M{z)\. In 
the limit of a strong anisotropy, the magnetic state could 
be viewed as parallel planes of magnetic solitons, where 
each soliton twists the magnetization- vector by 90°. 



We will now address the question, how superconduc- 
tivity nucleates in the presence of such a filamentary 
magnetic structure. From our analysis in the previous 
section, we know that in the limit I^^Llfco > v2 a ho- 
mogeneous superconducting order nucleates despite the 
modulated magnetic background. In the opposite case, 
IClI^o < v2, for which a spatially modulated supercon- 
ducting order is found, first of all the solitons in the 
magnetization lift the continuous degeneracy of the solu- 
tion ([9]) parameterized by ip. 

Then, the question arises whether the superconduct- 
ing order parameter nucleates at the soliton (domain 
wall) position or rather in the interior of the domains. 
To address this, let us consider the nucleation of bound 
states at an isolated domain wall at which the mag- 
netization is tilted from Tn{r) \\ b^ for z — >■ — oo to 
m{r) II By for z -^ oo. We shall for the moment assume 
that the amplitude |m(r)| = |M(z)| = toq is constant 
across the domain wall. The magnetization can thus be 
parametrized by a single function 6(z) = —9{—z) with 
9(z -^ ±oo) = ±7r/4 as 
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Insertion of the magnetization ( |14a[ ) in the free energy 
for the superconducting order parameter, Eq. ([3]), yields 
upon Fourier transformation in the x and y coordinates 
the variational equation 
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- sm[0{z) + tt/4] - cos[e{z) + 7r/4] . (14c) 



Here, we have chosen kx = ky = k in accordance with 
the symmetry of the problem and to obtain a binding po- 
tential. Eq. ( 14b ) is reminiscent of the one-dimensional 



Schrodinger equation with the potential V{z) and an en- 
ergy eigenvalue given by the square bracket on the rhs. 
This analogy immediately delivers the inequality 
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since the energy of the lowest bound state is always larger 
than the potential minimum V{Q) — ~\f2. Via the tem- 
perature dependence of a, the energy eigenvalue of the 
lowest bound state determines the temperature at which 
the superconductivity nucleates at the domain wall. The 

pdw 

satisfies 



transition temperature T^n of the bound state therefore 
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where T^n is the bulk transition temperature defined in 



sc 



eq. ( 4b ) . This inequality holds independent of the con- 



crete form of the function Q{z). We conclude that super- 
conductivity will nucleate in the interior of the domains 
rather than at the magnetic solitons (domain walls). 



This result relies on the assumption that |m(r)| is spa- 
tially constant. If we relax this assumption and consider 
the case in which the magnetization is suppressed near 
the domain wall to a value TOq < toq, the superconduct- 
ing order near the domain wall would be less affected by 
the paramagnetic depairing effect, represented by th e pa- 
rameter u. The upper bound for T^q in eq. (15b) then 



-(1) 



exceeds Tg^j , thereby opening the way for an reversion of 
the inequality T|j < T. 



.(1) 
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V. CONCLUSIONS 

In our study on the coexistence of ferromagnetism and 
superconductivity in noncentrosymmetric materials we 
found that ASOC, represented by Lifschitz-terms in the 
free energy expansions, gives rise to unusual modulations 
of the order parameters. Where a centrosymmetric ma- 
terial would have a homogeneously magnetized ferromag- 
netic ground state, an arbitrarily small parity violation 
generates a state of helical magnetization. 

The possible modulations in the superconducting or- 
der parameter in presence of magnetic order were found 
to be governed by the ratio of two length scales, |^l| and 
the characteristic length of magnetic modulations k^ . 
When the magnetic length scale is larger, the supercon- 
ducting order is not homogeneous. This extends to the 
limit of homogeneous magnetization k^ — > oo, where a 
complex phase modulation was found. In presence of the 



helical magnetic phase, a state of simultaneous complex 
phase and amplitude modulation develops, if the mag- 
netic wavelength is large enough. 

We closed our discussion with the consideration of a 
magnetic state of filamentary solitons, which is obtained 
for a strong anisotropy parameter in the magnetic free en- 
ergy. We found that this will force the order parameter 
of superconductivity to nucleate with the same filamen- 
tary structure, creating a stripe-like state with maxima 
in between two solitons. 

The complex phase winding of the order parameter of 
non-centrosymmetric superconductors that are exposed 
to a magnetic field was noticed theoretically a long time 
ago^*, but an experimental verification of this state is still 
lacking. The inhomogeneous superconducting states that 
we find in this work show both complex phase and ampli- 
tude modulations. The latter would entail an anisotropic 
resistivity drop at the criticality for current directions 
perpendicular and parallel to the wavevector of the mod- 
ulation and are therefore more accessible to an experi- 
mental observation. In fact, such anisotropic supercon- 
ducting transitions were reported for the antiferromag- 
netically ordered noncentrosymmetric CeRhSis in ref. [71 
where superconductivity and magnetism show and inter- 
esting interplay. 
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